Dynamically stabilized decoherence-free states in non-Markovian open fermionic 

systems 



O 



m 
in 

o 

(N 



Heng-Na Xiong, Wei-Min Zhang0 and Matisse Wei- Yuan Tu 
Department of Physics and Cenert for Quantum information Science, 
National Cheng Kung University, Tainan 70101, Taiwan 

Daniel Brauifl 

Laboratoire de Physique Theorique, Universite Paul Sabatier, 
Toulouse III, 118 route de Narbonne, F-31062 Toulouse, France 
(Dated: July 9, 2012) 

Decoherence-free subspaces (DFSs) provide a strategy for protecting the dynamics of an open sys- 
tem from decoherence induced by the system-environment interaction. So far, DFSs have been pri- 
marily studied in the framework of Markovian master equations. In this work, we study decoherence- 
free (DF) states in the general setting of a non-Markovian fermionic environment. We identify the 
DF states by diagonalizing the non-unitary evolution operator for a two-level fermionic system at- 
tached to an electron reservoir. By solving the exact master equation, we show that DF states can 
be stabilized dynamically. 



I. INTRODUCTION 

Quantum decoherence is a fundamental issue in open 
quantum systems and is also the most important obstacle 
opposing the realization of a large scale quantum com- 
puter. Unavoidable interactions with noisy environments 
typically render initially prepared pure states mixed very 
rapidly. Nonetheless, a subspace of Hilbert space can ex- 
ist where a system undergoes a unitary evolution irre- 
spective of the interaction with its environment. Such 
a decoherence-free subspace (DFS) can in principle pro- 
vide a theoretically perfect strategy to protect a system 
against quantum decoherence. The possibility of DFSs 
was already pointed out by Zurck 1] , who observed that 
if the interaction Hamiltonian of the system with the en- 
vironment has a degenerate eigenvalue, then superposi- 
tions of the corresponding eigenstates remain coherent. 
Generally, degenerate eigenvalues can notably arise if the 
coupling has a certain symmetry, and the resulting pro- 
tection against decoherence has been observed early in 
the context of rotational tunneling ^^^^ 
1990s, DFSs were rediscovered independently by several 
groups p^l - fisj . The theory of DFSs has already been 
extensively discussed in the context of the symmetries 
of the Hamiltonian [l^, E^l i semigroup dynamics in the 
language of quantum master equation |14| - [l9{ . and the 
operator sum representation based on Kraus operators 
[201. Meanwhile, the existence of DFSs has also been 
verified experimentally with polarization-entangled pho- 
tons [2l| - [23t . trapped ions 0, HI], nuclear spins using 
nuclear magnetic resonance techniques (26l - [29| . and neu- 
tron interferonietry (soj . These experiments show that 
encoding quantum information in a DFS can significantly 
prolong the storage time. Therefore, DFS have attracted 
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wide interest for applic ations in fault-tolerant quantum 
computatio n. llTl ISll . I33l438l | , long distance quantum com- 
munication |39| ■ quantum key distribution [40[ , quantum 
teleportation |4j|. quantum metrology [43 [ , robust quan- 
tum repeaters [43l. |44[ . and coherent quantum control 
[45l - l47j . The requirement of a symmetry in the cou- 
pling to the environment is not always necessary psl ). 
Recently, a new Heisenberg-limited metrology protocol 
was proposed which exploits the evolution of a DFS due 
to a collective change of the couplings to the environment 

According to Refs. [3, [l^, a DFS in a system ruled 
by Markovian dynamics is defined formally by the van- 
ishing of the nonunitary Lindblad (decohering) part in 
the Markovian master equation [53|, see Eq. This 
defines a Hilbert subspace in which dynamics is locally 
(in time) unitary. Even so, decoherence can still arise 
if the Hamiltonian of the system drives a state out of 
the DFS. This happens on the time scale of the system 
Hamiltonian that is typically much longer than the mi- 
croscopic decoherence times [s^l- However, leaking out 
of the DFS can be suppressed by coupling the system 
relatively strongly to the environment [5l|. An al- 
ternative definition of a DFS in the Markovian context 
has been given in These authors not only derived 

necessary and sufficient conditions for the vanishing of 
the nonunitary part of the master equation [see Eq. (j4])] , 
but also found criteria for globally DF states, for which 
the dynamics resulting from the full master equation (in- 
cluding the Hamiltonian part) remains unitary. 

The extension of such DFS criteria to non-Markovian 
dynamics is not straightforward. The general non- 
Markovian master equation of the Nakajima-Zwanzig 
form 52] involves a complicated time-non-local mem- 
ory integrand in the nonunitary terms. However, the 
exact master equations that describe the general non- 
Markovian dynamics have been recently developed for 
some class of open quantum systems, including quantum 
Brownian motion |53| - |55| . entangled cavities with vac- 
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uum fluctuations [56|, coupled harmonic oscillators jSTI . 

t quantum dot electronic systems in nanosturctures 
l60l | , various nano devices with time-dependent exter- 
nal control fields 'gT], nanocavity systems including ini- 
tial system-reservoir correlations [62|, and photonic net- 
works imbedded in photonic crystals [g^. These exact 
master equations can all have the nonunitary Lindblad 
form, but the decoherence rates are time-dependent and 
may become negative during the time evolution. This 
is diflFerent from the Markovian case, where the decoher- 
ence rates are always positive. Since the fact that all 
rates have the same sign plays an essential role in show- 
ing the necessity of the DPS criterion, other possibilities 
of creating a DF state may arise in non-Markovian cases. 

In this work, we study the dynamics of an open 
fermionic system by solving the exact non-Markovian 
fermionic master equation [59146 Ij . We find a new type 
of pure DF states which arise from the fact that certain 
time-dependent decoherence rates in the master equation 
switch themselves off after the system reaches a stable 
state. We call these DF states dynamically stabilized 
DF states. They are generated in particular for non- 
Markovian environments. The mechanism how the DF 
states arise is therefore very different from the known 
mechanism obtained from the Markovian master equa- 
tion. For a two-level fermionic system coupled to an 
electron reservoir, we find two dynamically stabilized DF 
states that possess full quantum coherence between the 
singly occupied states of the original two levels. Practi- 
cal applications of the dynamically stabilized DF states 
in electron spin and charge qubits for quantum informa- 
tion processing are expected. 

The paper is organized as follows. In Sec. II, we briefly 
review the general criterion of DFSs based on the Marko- 
vian master equation formalism. In Sec. Ill, we discuss 
the dynamics of electron systems in nanostructures via 
the exact master equation. In particular, we consider a 
two-level electron system whose two levels are coupled 
identically to an electron reservoir. We find that the 
system Hilbert space can be split into two closed sub- 
spaces. Then in Sec. IV, in terms of the full Lindblad 
generator, we discuss all the possibilities how DF states 
can arise in the system. We find that the vanishing of 
one of the two decoherence rates in the problem can give 
rise to dynamically stabilized DF states. With only one 
decoherence term in the master equation remaining, the 
well-known criterion of DFSs for Markovian decoherence 
still provides a necessary and sufficient condition for DF 
states. In Sec. V, we investigate under what conditions 
the dynamically stabilized DF states can be reached by 
the same dissipative process that allows their existence. 
We find that the initial state and the details of the dy- 
namics given by time-dependent non-equilibrium Green 
functions determines which DF state can be reached. In 
Sec. VI, we discuss the generation of DF states in the 
Born-Markovian (BM) dynamics within our exact frame- 
work. We show that the DF state in the BM dynamics is 
a special case of the exact solution. Finally, a conclusion 



is given in Sec. VI. 



II. CRITERION FOR DFS IN MARKOVIAN 
CASE 

In this section, we will briefiy review the general cri- 
terion of DFS based on the Markovian master equation 
following Ref. Consider an open system S coupled 

to an external noisy environment E. One can write the 
total Hamiltonian as 



H ^ Hs (E) Ie + Is He + Hj, 



(1) 



where / is the identity operator and Hj denotes the inter- 
action Hamiltonian between the system and its environ- 
ment. If the dynamics of the system-bath is Markovian 
and the system and the environment are initially decou- 
pled, the master equation for the density matrix of the 
system takes the Lindbladian form |65l - [67| 



[Hs,p{t)]+L[pit)] 



(2a) 



N 



^ W = I E «a {2F^pFl - FlF^p - pFiF^) , (2b) 



where Hs = Hs + A is the renormalized system Hamilto- 
nian with A a possible hermitian contribution from the 
environment ("Lamb shift"), F^ are orthogonal opera- 
tors on the system Hilbert space 'Hs, and denote real 
positive coefficients. Thus, the commutator involving Hs 
in Eq. ([2a]) determines an effectively unitary evolution of 
the system, while the decohering effect induced by the 
environment is totally accounted for by the non-unitary 
term L [p] . 

For the Markovian master equation ([5]), the condition 
of an instantaneous DFS at time t amounts to the van- 
ishing of L [p] , that is [11] , 



L [PDF {t)] = 0. 



(3) 



This ensures that pdf (i) obeys a unitary evolution 
PDF {t) — —i Fls, puF {t) at time t. This does not imply 
necessarily unitary evolution at all times, as the evolu- 
tion due to Hs can drive the system out of the DFS. In 
[ill the DFS found from © was therefore called DFS to 
the first order. 

A sufficient and necessary condition for a state |A:df) 
to be locally DF is given in terms of the operators {Fa} 

M by 



Falkop) = Cq, I fcDF),Va, fcDF 

N N 

and ^aaFlFalki^p) = X! ""'"^"'^'^df 



(4) 



Eq. ([U means that the DF states are degenerate eigen- 
states of all the operators {Fa}, and degenerate eigen- 
states of ^a=i o-aFaFl. A special case is given by = 
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Va, in which case the second condition in ^ is auto- 
matically fulfilled. In [l^ it was shown that the space 
spanned by the |A;df) is a DFS at all times, if and only 
if in addition to satisfying Q it is also invariant under 

Hs. 

The generalization of Q to non-Markovian master 
equations is not straight-forward. However, it was re- 
cently shown f7(J| that under certain conditions the solu- 
tion of a non-Markovian master equation [52j with a finite 
time memory kernel can be at the same time a solution 
of a local-in-time non-Markovian master equation. In 
principle, local in time generalizations of the Markovian 
master equation (j2bp may be obtained by making both 
the rates Oq. and the operators Fa time-dependent. Inter- 
estingly, the exact master equations describing the gen- 
eral non-Markovian dynamics for a large class of bosonic 
and fermionic systems [H, [59l - [63j have been developed 
recently. They all have the Lindbladian form of Eq. ([2|) , 
but the decoherence rates Oq, in the nonunitary term (|2bp 
are time-dependent and local in time, whereas the op- 
erators Fa are time- independent. The time-dependent 
decoherence rates are determined microscopically and 
nonperturbatively by the retarded and correlation Green 
functions in non-equilibrium Green function theory (63 [. 
where the back-actions from reservoirs are fully taken 
into account. As a result, the non-Markovian dynamics 
are fully characterized in terms of the time-non-local re- 
tarded integrand in the Dyson equatioii, which governs 
the nonequilibrium Green functions (68l l69j . In addition 
to being time-dependent, the decoherence rates can be- 
come negative for short times, representing in a certain 
sense the back-flow of information from the environment 
to the system 59, 71 74]. Since the fact that all have 
the same sign is an important requirement for showing 
that (m is necessary for a Markovian DF state (see the 
argument after eq.(3.15) in [l^), additional DF states 
may arise in the non-Markovian case. In the following, 
we will provide a new mechanism for generating dynami- 
cally stabilized DF states in fermionic systems, based on 
the exact fermionic master equation developed recently 

MM. 



III. EXACT MASTER EQUATION 

We consider a general nanoelectronic system with N 
energy levels coupled to an electron reservoir. The Hamil- 
tonian for the system, the electron reservoir, and the in- 
teraction between them read 



N 



Ck, 



i=l 
N 



i=l k 



Here a\ and at are electron creation and annihilation op- 
erators for ith level with energy e^. The operators cj. 



and Cfc denote electron creation and annihilation oper- 
ators for the energy level Sk of the electron reservoir. 
The coupling strength between the system and the reser- 
voir is described by the Vik G K, with the phase (t>i of 
the coupling made explicit. The total hamiltonian is 
H = Hs + Hb + Hj . 

The exact master equation for such system was ob- 
tained in 



p{t) = -i Hs{t),p{t) +L[p{t)] 



(6a) 



where the renormalized system Hamiltonian Hs (t) and 
the decoherence term L [p (t)] take the form 



N 



(6b) 



N 



L[pit)]= {'^yW[2ajP(i)al-a|ajp(t)-p(i)ajaj] 

+K^j (t) [2alp{t)aj ~ ajajpit) - p{t)ajal] } . (6c) 

Here the shifted eij{t) and decoherence rates Kij{t) and 
Kij(t) are all time-dependent but local in time. They 
are determined microscopically and nonperturbatively in 
terms of the retarded and correlation Green functions by 
eliminating completely all the reservoir degrees of free- 
dom (i.e. tracing over all the states of the environment). 
Their explicit forms are shown in Eq. (9) of Ref. [59| . 

To be more specific, let us consider the case of a nano- 
electronic system with N = 2. Physically, such a system 
may be realized by a double quantum dot system in which 
each dot has a single active energy (on-site) level, cou- 
pled to electrodes with all spins polarized in both the dots 
and the electrodes. Another example for = 2 is given 
by a single-level quantum dot coupled to electrodes with 
allowed spin flips between two antiparallel directions. 

Furthermore, we assume that the two energy levels of 
the system are degenerate: ei = €2 = cq. Practically, the 
energy degeneracy is easier to be realized in the second 
setting than in the first. We also assume that both levels 
have the same coupling strength to the electron reservoir, 
i.e., Vik — V2k — Vfe/%/2. Then we introduce two effective 
fermion operators 

A+ = (ai + e"*a2) , A- ^ (-e~''^ai -I- 02) , 
V 2 V 2 

with (j) — (pi ^ (t>2- The Hamiltonian ([5]) can be rewritten 
in terms of A± as follows 

Hs ^ e+AlA+ + e^AlA^, Hb =^£fc4cfc, 

k 

k 

where the effective energy levels e± = eo are still degen- 
erate. 
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The system Hamiltonian is diagonalized in terms of 
A±, such that the original system is equivalent to an 
effective system which has two decoupled energy levels 
e±, out of which only one energy level (e-|_) couples to 
the electron reservoir. As a result, the corresponding 
exact master equation becomes 



|±) = ^±|v) are superpositions of singly occupied states 
of the original two coupled levels, 



p{t) = -i Hsit),p{t) +ib(t)]. 



(8a) 



with the new expressions of Hs (t) and L [p (t)] 

Hsit) ^e+it)AlA+ + e_AlA_, (8b) 
L [p (t)] = ^ (t) [2A+P (t) Al - AlA+p (t) - p (t) AlA+ 

+ K (t) \2AIp (t) A+ - A+AIp (t) - p {t) A+AI 



(8c) 



The decoherence rates are K{t) — ^ (t) — 2^ and k {t) = 

^j^. The renormalized energy level £+ (t) and the rates 
7 (t) and 7 (t) are determined exactly by 

I+{t) = -lm[u{t)u-\t)], 



7(t) = v{t) - 2v{t)Re[u{t)u-^t)], 



(9a) 
(9b) 
(9c) 



and u{t) and v (t) are the retarded and correlation 
Green functions in the Schwinger-Keldysh nonequilib- 
rium Green function theory [Glj . They obey the integral- 
differential equations 



dt 

vit) 



d r* 

—u {t) + ie+u {t)+ g{t-T)u (r) dr = 0, (10a) 

J to 

rt nt 

dTi / dT2U (ri) g (t2 - n) u* (T2) , (10b) 

fto Jto 



subject to the boundary conditions u (to) = 1 . The 
integration kernels read 5 (t) = i^) ^"^'^^ 

and g{T) = J (w) / (w) e"*'^^, where the spec- 

tral density of the reservoir is given by J (w) = 
27r \Vk f 6 (uj-uJk), e^nd fiu) = 1/ (e^^""^) + l) is 
the initial electron distribution of the reservoir. 

Since the effective energy level e_ is apparently decou- 
pled from the electron reservoir, one may naively think 
that the state generated by the operator A^_, namely 

A^jO) becomes naturally a DF state. This is actually 
not true. The term with K{t) in ([5c)) will drive this 
state into another state since A^^_A^1|0) ^ 0. However, 
there is an occupation constant of motion in this system, 

A''_A-,H = 0. This symmetry separates the system 

Hilbert space into two closed subspaces 71+ = {|v), |-|-)} 
and H- = {|— ),|d)}, corresponding to the occupation 
7V_ = (ALA_) = and 1, respectively. Here |v) (|d)) 
is the vacuum (doubly occupied) electron state, while 



(|l>+e- 
(-e'^ll) 



|2))/V2, 
|2)) /V2, 



(11) 



and \i) = aj|v) (i = 1,2). The two states |±) are the 
superpositions of the two original energy levels with the 
relative phases (j) and (j) + t:, where (j) is arbitrary. As a 
consequence, starting from any initial state in the closed 
subspace ("H-), the system will be kept in this sub- 
space throughout the evolution process. 

Eq.® has the standard form of Lindblad master equa- 
tion, except that the decoherence rates, K(t) and Ti{t), 
can depend on time and even become negative. They are 
local-in-time and determined microscopically and non- 
perturbatively from Eq. (ITO|) . The DFS criterion of 
Eq. (U) is then still sufficient, but may not be necessary. 
The operators A\_ and ^+ act on the four basis states 
{|v), |— ), |-|-), |d)} according to the following relations: 





|v) = |+), 




|v)=0. 




l+) = o, 


A+ 


1+) = |v 


Al 


|-> = |d>, 


A+ 


|->-o, 


At 


Id) =0, 


A+ 


Id) = 1- 



One checks that both A]_ and A^ have a four-fold 
degenerate eigenvalue Cq = 0, but the corresponding 
eigenspaces are only two-dimensional, given by {|+), |d)} 
for A'\_ and {|— ), |v)} for A^. Since these two spaces do 
not overlap, condition @ is not satisfied, and as long as 
both K{t) and K{t) are positive, there is no DFS. In the 
next section, we explore whether additional DF states 
can exist if K{t) and K{t) are not both positive, and show 
that a DFS can arise dynamically, when at least one of 
the two rates K{t) and K{t) vanishes. 



IV. DF STATES 

A. Eigenstates of the Lindblad operator 

Since, as discussed above, the Markovian DFS crite- 
rion ^ may not be necessary for a non-Markovian sys- 
tem, we define a local (in time) DFS by the vanishing 
of the non- unitary term L[p{t)]. Clearly, this leads to 
local unitary time-evolution. In the following, we will 
discuss all possibilities how L[p{t)] — can arise by cal- 
culating the eigenvalues of the full Lindblad generator L 
in Liouville space. If L has a zero eigenvalue, the corre- 
sponding eigenstate is DF. Conversely, a local DF state 
is by definition an eigenstate of L with eigenvalue zero. 
However, as the eigenstates of L do not necessarily have 
all the properties of a density matrix, such states may 
not be physical. One must therefore examine for each 
eigenstate whether it is a physical state or can be com- 
bined with other eigenstates corresponding to the same 
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(degenerate) eigenvalue to form a physical state. In this 
way one can find all physical DF states possible. 

In the basis of {|v), |+), |— ), |d)}, the density matrix 
takes the general form 

p{t) =p,,{t) |v)(v|+p++ 

+ p+-{t)\+){-\ + pX- Wl-X+I 

+ (i)|->(-|+Pdd(<)|d)(d|. (12) 

Coherences between states with different particle num- 
bers are not permitted due to the particle number 



super-selection rule. Since there are only six nonzero 
elements of p(i), we can define the basis {|1) = 
|v)(v|, |2) ^ |3) ^ |+)(-|, |4) ^ |-)(+|, 

|5) = |— )(— I, |6) = |d)(d|} in Liouville space, 
which is orthogonal with respect to the scalar prod- 
uct {A\B) = trA^B. In this basis, the density ma- 
trix can be rewritten as a column vector, \p (t)) = 

[pvv (t) , P++ (t) , p+_ it) , p_+ {t) , p__ {t) , pdd {t)f. The 
Lindblad operator L is represented by a matrix Lf, and 
its action on a density matrix reduces to a simple matrix 
multiplication of Lt with \p(t)), where Lt reads 



Lt 



-2k {t) 
2k (t) 







2k {t) 
-2k (t) 














{K{t) + K{t)) 












(k it)+Kit)) 













-2k (t) 
2k (t) 








2k (t) 
-2k (t) 



(13) 



The eigenvalues and the corresponding eigenstates of Lt 
are 



h = -2 [k {t) + K{t)] 

^3 = - [« {t) + K {t)] , 
h^-[K (t) + K {t)] , 
^5 = -2[k {t)+K{t)] 

k = o, 



and 



(14) 



h) 
h) 

h) 
k) 



K{t) + K (t) 

^[-|1> + |2)], 



[Kit)\l)+K{t)\2)], 



= |3), 
= |4), 
1 . 



-|5) + |6)], 



K{t) + K (t) 



[Kit)\5)+Kit)\6)]. (15) 



We see that all eigenvalues come in pairs, and there are 
always at least two eigenvalues equal to zero. In writing 
([15]) we have assumed that K{t) + K{t) ^ 0. The case 
of K(t) 4- K(i) = will be discussed below. The normal- 
ization used for |Zi) and Zg) is convenient, as in this way 
these two states can be interpreted directly as density 
matrices if both K{t) and 'kit) are positive. If one of the 
rates is negative (and the other non-zero), both states 
become non-positive and therefore cease to be physical 
states. Moreover, since they have orthogonal support. 



no linear combination of them can bring about a positive 
state. Similarly, states [^2) and IZ5), as well as any linear 
combination of them, are clearly non-positive. Finally, 
states 1^3) and IZ4) as well as any linear combination of 
them are traceless and are therefore not physical states 
either. Since I/2), . . . , IZ5) are independent of K(i) and K(t) 
(and are therefore never physical states, regardless of the 
values of K{t) and K(t)), the only possibility of having a 
physical DF state is through |Zi) and with both k[C) 
and K(t) non-negative, or one of them vanishing (in the 
latter case one may always choose the eigenvector with 
positive global sign). 

If K (t) = and K(t) ^ 0, and \l^ are two pure 
DF states |2) and |6), that is, the subset of {|+), |d)} 
contains all the possible DF states in this case. Likewise, 
if K (<) = and k (t) 7^ 0, and are the two pure DF 
states |1) and |5), i.e., {|v),|-|-)}. Below, by examining 
an explicit example, we will show that the vanishing of 
one of the decoherence rates is physically feasible after 
some time ts, when the system reaches its steady state, 
see Fig. ^ 

If K{t) > and K{t) > 0, the two eigenstates \li) 
and are time-dependent mixed states. They are 
decoherence free as much as they are locally stationary 
states due to local detailed balance. From a perspective 
of application for quantum information processing, these 
states are less interesting. They are the analogues of 
thermal equilibrium states that are stationary under a 
Markovian relaxation process. 

It remains to consider the case K{t) + K{t) = 0. All 
eigenvalues vanish, but as long as K{t) 7^ the eigenvec- 
tors are still given by Eq. with the only difference 
that the normalization of \li) and j/g) through the pref- 
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actor 1/(k(<) + K{t)) has to be removed. \li) becomes 
cohnear with [^2), and \Iq) cohnear with I/5). The dimen- 
sion of the space of eigenvectors of Lt is reduced to four, 
and Lt can therefore not be fully diagonalized. As the 
linearly independent eigenvectors \li), . . . , IZ4) are never 
physical states, this means that Lt has no eigenstates 
that are physically possible, and therefore no DF states 
exist. 

If K{t) = K{t) = at some time (for example, if the re- 
tarded Green function u (t) took a nonzero steady value, 
it would be possible that both rates vanish for i > tg), we 
have Lt = 0, and the whole Hilbert space becomes a dy- 
namically stabilized DFS. From a quantum information 
perspective this would be of course an ideal situation. 
Unfortunately, in the fermionic system interacting with 
an electron reservoir considered here, it appears that this 
situation does not arise, as shown in Fig. [51 

In summary, as long as we restrict ourselves to pure 
states for our non-Markovian master equation and dis- 
card the "trivial" case Lt = 0, the DFS is still given en- 
tirely by the Markovian criterion as pure DF states 
only exist if exactly one of the two terms (proportional 
to either hi{t) or Hit)) in the Lindblad superoperator re- 
mains, and the logic of the proof of necessity of condition 
remains intact in such a situation. However, the time 
dependence of K{t) and k(<) brings about a new freedom, 
and allows for the dynamical stabilization of DF states 
through the switching off of one of the decoherence rates. 



B. Physical realization of DF states 

In the following, we will discuss to what extent the DF 
states just discussed can be reached through the same 
non-Markovian dynamics described by Lt- We consider 
the double quantum dot system from above coupled to an 
electron reservoir with a spectral density of the Lorentz 
form [silTllli, 



J{uj) = 



(16) 



where F is the system-reservoir coupling strength, and d 
the bandwidth of the effective reservoir spectrum. In the 
well-known wide-band limit, i.e., d — !■ 00, the spectral 
density approximately becomes a constant one, J{uj) — >■ 
F. This corresponds to the Markovian limit. 

With the above spectral density, the solution of u{t) 
obeying Eq. (|10al) can be obtained analytically 



u{t) 



2dr 
[1- 



dp e 



-dr 



dtl 



=-(«++!)* 



,d^2F, 
d = 2F, 



(17) 



where dr — \/(P~2rd and dp = d ± dr- Obviously, 
after some time tg, u (t) always decays to zero, as shown 
in Fig. [T] where the different behaviors of the amplitude 
of u (t) corresponds to different bandwidths d. For the 




FIG. 1: The exact solution of \u{t)\ for an electron reservoir 
with the Lorentz spectral density. Here we take eo = 0.2r. 
Note that \u{t)\ decays exponentially for large bandwidth d, 
corresponding to the weakly non-Markovian dynamics. For a 
small bandwidth d < 2r, the strong non-Markovian memory 
effect brings the short-time oscillations for |?t(f)|. 



bandwidth d > 2r (weakly non-Markovian case), u{t) 
exponentially decays to zero, which is a result similar 
to the Markovian dynamics, see the discussion in Sec. 
VI. When the bandwidth d < 2r, the non-Markovian 
memory effect of the reservoir becomes significant, which 
induces a short-time oscillation for \u {t)\. 

Using the solution of Eq. (IT7)) . combined with 
Eqs. (|10b[l and ([9]), we can easily calculate the decoher- 
ence rates k (t) and k (t) in eq. The result is plotted 
in Fig. [5] for the cases of weakly (d = lOF) and strongly 
(d ~ 0.5F) non-Markovian dynamics. In experiments, 
one can adjust the external bias (^) to raise or lower the 
Fermi surface of the electron reservoir. Here we display 
three cases where the bias is much higher (/z = lOF), rel- 
atively small (p — 0), and much lower (/x = — lOF) than 
the quantum dot energy level eg. First, we see that when 
d = 0.5F, as shown in Fig. [5] (c) and (d), k (t) and k (t) 
can jump from a positive value to a negative value repeat- 
edly during the evolution process, which corresponds to 
the forth- and back-flows of the information between the 
system and the environment, in evidence of the strong 
non-Markovian dynamics. 

We observe that, no matter what the values of the 
width d and the bias /j, are, the two decoherence rates 
K{t) and K{t) never satisfy K{t) + /5(t) — 0. However, 
if we apply a large bias to raise (or lower) the Fermi 
surface of the electron reservoir much higher (or much 
lower) than the dot level, one of the two decoherence 
rates is switched off after a time scale of a few 1 /F, which 
implies the existence of a DFS. 

As shown in Fig. [2] (a) and (b), in the weakly non- 
Markovian case, when applying a positive bias = lOF, 
the decoherence rate k (t) shows a positive peak in the 
beginning and then decays rapidly to a zero steady value 
on a time scale of a few 1/F, while K{t) turns negative 
first, and then climbs to a nonzero steady value. For 
the strongly non-Markovian regime (d = 0.5F), the same 
situation happens, see the red curves in Fig. [21(c) and (d). 
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(a) d=lor 



(b) d=lor 



tively given by 




FIG. 2: The exact solutions of decoherence rates K{t) and 
K{t) for different external bias voltage fi = eV. (a)-(b) for the 
weakly non-Markovian dynamics (with d = lOF); (c)-(d) for 
the strong non-Markovian dynamics (with d — O.SF). Here 
we set eo = 0.2r, ksT = 0.3F. We note that generally, K{t) 
and K{t) are non-zero during the time evolution. However, 
for a large positive (or negative) bias voltage, one of them is 
switched off after some specific time. 



The decoherence rate reaches K{t) = on a time scale 
ts of a few l/F, while /5(<) keeps jumping from positive 
values to negative values repetitively. In this case, the 
dynamics for t > tg is described by the following master 
equation 



Hs (t) = e+ (t) AU+ + e-^U-, 



(18a) 
(18b) 

L [p (t)] - H (t) [2AIp (t) A+ - A+AIp (t) - p (t) A+Al 

(IBcj 



Then the states |-f ) and |d) become possible dynamically 
stabilized DP states, since AI\+) = and also ^^|d) = 
leads to the vanishing of the nonunitary part L[p{t)] after 
t > ts- Therefore, the Markovian DFS criterion of Eq. (|4]) 
still is both necessary and sufficient for a time-local DFS 
when only one decoherence rate is turned on. 

On the other hand, if we apply a negative bias to 
the electron reservoir, e.g. p = — lOF, as shown by the 
black curves in Fig. [21 we find that for both weakly and 
strongly non-Markovian dynamics, the decoherence rate 
K(t) goes to zero very quickly, while K{t) either reaches 
a nonzero steady value (in the weakly non-Markovian 
regime), or keeps jumping from a positive value to a neg- 
ative value for all times (in the strongly non-Markovian 
regime). Then the master equation for t > is effec- 



p it) = -I ^^Hs (t) , p (t)J +L[p (t)] , (19a) 

Hs{t) ^e+{t)AlA++e^AlA^, (19b) 

L [p (<)] = (t) [2A+P {t) 4 - AlA+p (<) - p it) AlA+ 

(19c) 



Again, the DFS criterion Eq. is then both neces- 
sary and sufficient for DF states. Since ) = and 
A+jv) — 0, then we obtain that |— ) and |v) are possible 
dynamically stabilized DF states. 

In summary, the above results show that the vanishing 
of one of the decoherence rates implies that |d)} of 

{|— ), |v)} become dynamically stabilized DF states. 



V. PHYSICAL REALIZATION FOR 
DYNAMICALLY STABILIZED DF STATES 

In the following, by examining the general exact solu- 
tion of master Eq. ([8|), we will prove that all pure dy- 
namically stabilized DF states {|+),|d)} and {|— ), |v)} 
predicted above are physically realizable through the de- 
coherence process given by the same master equation. 
We also give the exact conditions for generating these 
DF states in terms of the initial state, and a condition 
on the Green's functions u{t) and v{t). 

First, by solving ([8|), we can exactly give the elements 
of p (t) in Eq. (fT2|) in terms of the initial p (to) and the 
functions of u (t) and v (t) as 

Pvv (i) = [l-v (t)]pvv (io) + [l-v{t)~ \u {t)\^]p++ {to) , 
P++ (t) V (t) p^ (to) + [v (t) + \u {t)f]p++ (to) , 
p+_ it) = u{t) e"-V+-(io), 

P— (t) = [1-V (t)]/,__ (io) + [1-Vit)- \U it)f]pM (to) : 



Pdd (t) = V (t) (to) + [v (t) + \u (t)|>dd (to) ■ 



(20) 



This expression for p (t) is valid for an arbitrary spec- 
tral density of the reservoir. This solution confirms the 
fact that 'H+ = {|v), \+)} and V.^ = {]-), |d)} are two 
independent closed subspaces. For any initial state in 
the subspace (^-), the system will be dynamically 
stabilized in this subspace. 

The general solution ([20| shows that if v (t) ~ 1 is sat- 
isfied when t > ts, the initial vacuum state |v) converges 
to the stabilized state |-|-), while ii v {t) + \u {t)\ =1 
is reached when t > ts, the initial singly occupied state 
|-|-) converges to the dynamically stabilized DF state |+). 
That is, the same dynamically stabilized DF state |-|-) can 
be generated from different initial states in the same sub- 
space of under different stabilization conditions for 
u{t) and v{t). Similarly, one can obtain the remaining 
dynamically stabilized DF states |v), |— ) and |d) from 
different initial states under different stabilization condi- 
tion, as shown in Table HI 
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TABLE I: The dynamically stabilized DP states for different 
choices of the initial state of the system, plus the different 
stabilization conditions. 



(a) cl=10r 



initial state 


stabilization condition 


DP state 


|v) 


v{t) = (or 1) 


|v) (or 1+)) 


1+) 


v{t) + \u{t)f = (or 1) 


|v) (or 1+)) 


|-> 


v{t) = (or 1) 


|-> (or |d)) 


Id) 


v{t) + \u{t)f = (or 1) 


h> (or Id)) 



Besides the initial states listed in Table HI another more 
general initial pure state is a superposition of |+) and |— ) , 
namely |$) = a|+) + /3|-) with |ap + |/3p = 1. In this 
case, as one can easily check from Eq. (pHl) . the resulting 
state becomes 



p,,{t) = [l-v{t)-\u{t)\^]\a\'', 
P++ (t) = [v{t) + \u{t)\']\a\\ 
p+-{t)=u{t)e''-' a/3*, 

Pdd {t)^v{t)\P\''. 



(21) 



In the case that both a and /3 are nonzero, one has to 
have v{t) — and |w(t)p = 1 to generate a stabilized 
(pure) DF state |$). This is impossible unless the sys- 
tem is totally decoupled from the environment from the 
beginning. In other words, except for a stabilized mixed 
state, no stabilized DF state can be obtained. As a con- 
clusion. Table U lists all the possible pure stabilized DF 
states in this system. The present results confirm the 
statement in the last section that the only possible pure 
dynamically stabilized DF states are |±) and \d). 

For the Lorentz spectral density, we display the func- 
tions oiv{t) and v{t) + \u{t)\'^ in Fig. [3] Interestingly, by 
comparing Fig. |3] (a)-(b) and (c)-(d), we find that in the 
strongly non-Markovian case {d = O.SF), applying a rel- 
atively small bias voltage (for example \p\ = 2T) makes 
the steady values of v [t) and v [t) + \u {t)\ quickly ap- 
proach 1 or 0, whereas this is not the case in the weakly 
non-Markovian regime [d = lOF), where v{t) — > 0.9 for 
/i = 2r (and v{t) 0.1 for p. = -2V). This indi- 
cates that the back-flow of information from the reservoir, 
due to a strongly non-Markovian memory of the reser- 
voir, helps the stabilization of the system in the states 
of {|-|-),|d)} or {|v),|— )}. This facilitates the physical 
realization of the dynamically stabilized DF states. 

Furthermore, we can now prove that under the stabi- 
lization conditions listed in Table U one of the decoher- 
ence rates vanishes when t > ts- In fact, K{t) and /5(t) in 




(c) cl=0.5r 






5 10 


15 




n 




1 X — 


(d) d=0.5r 
> 






PIG. 3: The exact solution of v{t) and v{t) + |w(f)|^ for an 
electron reservoir with varying the external bias voltage p — 
eV. (a)-(b) for the weakly non-Markovian dynamics (with 
d = lOr); (c)-(d) for the strongly non-Markovian dynamics 
(with d = 0.5r ). Here we set eo = 0.2r, KbT = O.SF. By 
applying a large positive (or negative) bias voltage, the steady 
value of v{t) and v{t) + \u{t)\'^ will approach 1 (or 0). 



Eq. can be further simplified as 



Kit) 



iu(t)r d i-v{t) 



2 dt \u{t)\ 
,{t)\^ d vjt) 
~2~dt\u{t)\^ 



2 ' 



(22a) 
(22b) 



It clearly shows that the condition v{t) = 1 (or {v{t) + 
|u(t)P = 1) implies K{t) = 0, and the condition v(t) = 
(or v{t) + |u(i)P = 0) indicates K{t) = 0, namely, one 
decoherence rate is turned off for t > tg. In this case, the 
Markovian DFS criterion still provides a necessary and 
sufficient condition for generating dynamically stabilized 
DF states. However, for the non-Markovian dynamics, 
we have shown in TableUthat which final state is realized 
depends both on the initial state and the details of the 
dynamics given by the Green functions u{t) and v{t). 



VI. COMPARISON TO THE BORN-MARKOV 
DYNAMICS 

In the previous Sections, using the exact master equa- 
tion, we have shown the generation of DF states for a 
system coupled to a non-Markovian reservoir. For com- 
parison, in the following, we discuss the corresponding 
results in terms of the BM master equation. 

The BM dynamics usually corresponds to the case 
where the coupling strength between the system and the 
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electron reservoir is very weak, and the characteristic 
correlation time of the electron reservoir is sufficiently 
shorter than that of the system. In such a case, the elec- 
tron reservoir has no memory effect on the evolution of 
the system. Then the solution of u{t), v {t) are reduced 

to (Ha-iil 

^;BM(t)- [l-e-^('+)*]/(6+), (23) 

where — e+ + {6e+)^^ with the energy shift 
((5e+)„M = V f^°° ^li^^ Note that for the wide-band 
limit, we simply have e+ = e+ and J{e+) = F. Sub- 
stituting Eq. p3| into Eq. (jO]), we obtain the constant 
rates 

«(t)-ij(6+) [1-/(6+)], 

«(t)-ij(6+)/(e+). (24) 

Here / (e+) is the fermion distribution function of the 
electron reservoir at the frequency e+, i.e.. 

Thus the exact master equation is reduced to the BM 
master equation, where the decoherence rates are time- 
independent. This gives the standard Lindblad form for 
the Markovian dynamics. Based on the general DFS cri- 
terion of Eq. Q , we see that there is in general no DFS 
for this system in the BM limit. 

However, if we apply a large positive bias fi = eV such 
that, (/i — e+) /ksT ^ 1, then / (e+) 1, which leads to 
K — > 0. Hs leaves the states |+) and |d) invariant during 
the time evolution. The relation A+|-|-) = (or A+|d) = 
0) guarantees L [p {t)] = for these states. The states 
{!+), |d)} are therefore DF states in the BM limit under 
large positive bias. Likewise, applying a negative large 
bias to the electron reservoir such that (e+ — /i) /KgT ^ 
1, then / (e+) — and K(i) — ^ 0. In this case, the states 
{|v), |-)} are the DF states in the BM limit. 

In conclusion, the states {|+), |d)} and {|v),|— )} 
are also DF in the BM limit if one of the decoherence 
rates is switched off by properly tuning the bias voltage 
on the electron reservoir. This result is consistent 
with the result in the non-Markovian case discussed 
above. The apparent difference is that the DF states 
in the BM limit seem to exist without the dynamical 
stabilization processes. However, this difference is not 
crucial in reality. As is well-known [65, 66], the BM 
master equation with the constant decoherence rates. 



Eq. (l24l) . is derived under the condition t -ry where 
Tr is the characteristic time of the reservoir [65|. In 
other words, a stabilization time scale tg has implicitly 
been used in deriving the BM master equation, such 
that the decoherence rates become time-independent for 
t > tg ^ Tr ■ Therefore, the concept of the dynamically 
stabilized DF states, based on the exact non-Markovian 
master equation, gives the generalized picture of DFS. 
It contains the Markovian DFS as a special case. 



VII. CONCLUSION 

In summary, we have investigated the DFS of a non- 
Markovian fermionic open systems, based on an exact 
non-Markovian master equation developed recently. The 
master equation has a nonunitary term of the standard 
Lindblad form, but the corresponding decoherence rates 
are time-dependent and local in time. They are deter- 
mined microscopically and nonperturbatively from the 
Schwinger-Keldysh nonequilibrium Green functions and 
fully account for the non-Markovian memory effect. 

As concrete example we have studied a fermionic sys- 
tem with two degenerate energy levels coupled identically 
to a fermionic reservoir. We find that the whole Hilbert 
space is split into two closed subspaces. For any initial 
state in one of the subspaces, the system will remain in 
this subspace forever. By diagonalizing the full Lindblad 
operator we found that physical DF states exist if and 
only if one of the two relevant decoherence rates switches 
itself off dynamically. Such a situation can be achieved 
physically. Two of the DF states are coherent superpo- 
sitions with an arbitrary relative phase between the two 
original energy levels, which may be of physical interest 
for quantum computation. 

Which DF state is reached as result of the dissipa- 
tive dynamics depends both on the initial state and 
the details of the dynamics, as expressed by the time- 
dependent non-equilibrium Green's function. We show 
this explicitly by solving exactly the non-Markovian mas- 
ter equation. Interestingly, we find that the strongly non- 
Markovian memory can help to stabilize the DF states 
compared to the Markovian case. 



Acknowledgments 

This work is supported by the National Science Coun- 
cil of ROC under Contract No. NSC-99-2112-M-006-008- 
MY3, the National Center for Theoretical Science of Tai- 
wan. 



[1] W. H. Zurek, Phys. Rev. D 26, 1862, (1982). 



[2] E. B. Wilson, J. Chem. Phys, 3,276 (1935). 



10 



[3] H. C. Higgins, Mol. Phys. 7, 445 (1963). [35 
[4] J. H. Freed, J. Chem. Phys. 43, 1710 (1965). 
[5] K. H. Stevens, J. Phys. C 16, 5765 (1983). [36 
[6] W. Hiiusler and A. Hiiller, Z. Phys. B 59, 177 (1985). 
[7] A. Wiirger, Z. Phys. B 76, 65 (1989). [37 
[8] D. Braun and U. Weiss, Z. Phys. B 92, 507 (1993). 
[9] D. Braun and U. Weiss, Physica B 202, 264, (1994). [38; 
[10] G. M. Palma, K. A. Suominen and A. K. Ekert, Proc. 

Roy. Soc. London Ser. A, 452, 567 (1996). 
[11] L. M. Duan and G. G. Guo, Phys. Rev. Lett. 79, 1953, [39 
(1997). [40; 
[12] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 

(1997). [41 
[13] P. Zanardi and M. Rasetti, Mod. Phys. Lett. B 11, 1085 

(1997). [42 
[14] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 

Lett. 81, 2594 (1998). [43 
[15] D. Braun, P. A. Braun, and F. Haake, Lecture Notes in 

Physics 538, 55 (2000) (Proceedings of the 1998 Bielefeld [44 
Conference on Dccoherencc: Theoretical, Experimental, 
and Conceptual Problems"). [45 
[16] P. Zanardi, Phys. Rev. A 57, 3276 (1998). 
[17] D. Bacon, D. A. Lidar, and K. B. Whaley, Phys. Rev. A [46 
60, 1944 (1999). 

[18] R. Karasik, K. P. MarzUn, B. C. Sanders, and K. B. [47 

Whaley, Phys. Rev. A 77, 052301 (2008). 
[19] A. Shabani and D. A. Lidar, Phys. Rev. A 72, 042303 [48 

(2005). 

[20] D. A. Lidar, D. Bacou, and K. B. Whaley, Phys. Rev. 

Lett. 82, 4556 (1999). [49; 
[21] P. G. Kwiat, A. J. Berglund, J. B. Altepeter, and A. G. [50^ 

White, Science 290, 498 (2000). 
[22] J. B. Altepeter, P.G. Hadley, S.M.Wendelken, A. J. 

Berglund, and P.G. Kwiat, Phys. Rev. Lett. 92, 147901 [51 

(2004). 

[23] T. Yamamoto, K. Hayashi, A. K. Ozdemir, M. Koashi, [52 
and N. Imoto, Nature Photon. 2, 488 (2008). 

[24] D. Kielpinski, V. Meyer, M. A. Rowe, C. A. Sackett, W. [53 
M. Itano, C. Monroe, and D. J. Wineland, Science 291, 
1013 (2001). [54 

[25] C. Langer, R. Ozeri, J. D. Jost, J. Chiaverini, B. De- [55 
Marco, A. Ben-Kish, R. B. Blakestad, J. Britton, D. B. [56 
Hume, W. M. Itano, D. Leibfried, R. Reichle, T. Rosen- 
band, T. Schaetz, P. Q. Schmidt, and D. J. Wineland, 
Phys. Rev. Lett. 95, 060502 (2005). [57 

[26] L. Viola, E. M. Fortunato, M.A. Pravia, E. Knill, R. 

Laflamme, D.G. Cory, Science 293, 2059 (2001). [58 

[27] E. M. Fortunato, L. Viola, J. Hodges, G. Teklemariam, 

and D.G. Cory, New J. Phys. 4, 5 (2002). [59 

[28] J. E. OUerenshaw, D. A. Lidar, and L. E. Kay, Phys. 

Rev. Lett. 91, 217904 (2003). [60 

[29] D. Wei, J. Luo, X. Sun, X. Zcng, M. Zlian, and M. Liu, 

Phys. Rev. Lett. 95, 020501 (2005). [61 

[30] D. A. Pushin, M. G. Huber, M. Arif, and D. G. Cory, 

Phys. Rev. Lett. 107, 150401 (2011). [62 

[31] D. Bacon, J. Kempe, D. A. Lidar, and K. B. Whaley, 

Phys. Rev. Lett. 85, 1758 (2000). [63 

[32] W. T. Strunz, F. Haake, and D. Braun, Phys. Rev. A 67, [64 
022101 (2003). 

[33] A. Beige, D. Braun, B. Tregenna, and P. L. Knight, Phys. [65 

Rev. Lett. 85, 1762 (2000). 
[34] A. Beige, D. Braun, and P. L. Knight, New J. Phys. 2, 

22 (2000). [66 



Z. W. Zhou, B. Yu, X. Zhou, M. J. Feldman, and G. C. 
Guo, Phys. Rev. Lett. 93, 010501 (2004). 
L. A. Wu, P. Zanardi, and D. A. Lidar, Phys. Rev. Lett. 
95, 130501 (2005). 

P. Xue and Y. F. Xiao, Phys. Rev. Lett. 97, 140501 

(2006) . 

T. Monz, K. Kim, A. S. Villar, P. Schindler, M. Chwalla, 

M. Riebe, C. F. Roos, H. Haffner, W. Hansel, Phys. Rev. 

Lett. 103, 200503 (2009). 

P. Xue, Phys. Lett. A 372, 6859 (2008). 

Z. Q. Yin, Y. B. Zhao, Z. W. Zhou, Z. F. Han, and G. 

C. Guo, Phys. Rev. A 77, 062326 (2008). 

H. Wei, Z. J. Deng, X. L. Zhang, and M. Feng, Phys. 

Rev. A 76, 054304 (2007). 

C. F. Roos, M. Chwalla, K. Kim, M. Riebe and R. Blatt, 
Nature 443, 316 (2006). 

A. Klein, U. Dorner, C. M. Alves, and D. Jaksch, Phys. 
Rev. A 73, 012332 (2006). 

U. Dorner, A. Klein, and D. Jaksch, Quant. Inf. Comp. 

8, 0468 (2008). 

P. Cappellaro, J. S. Hodges, T. F. Havel, and D. G. Cory, 
Laser Phys. 17, 545 (2007). 

M. Kiffner, J. Evers, and C. H. Keite, Phys. Rev. A 75, 
032313 (2007). 

W. Wang, L. C. Wang, and X. X. Yi, Phys. Rev. A 82, 
034308 (2010). 

D. Braun, "Dissipative Quantum Chaos and Decoher- 
ence" , Springer Tracts in Modern Physics 172, Springer, 
Heidelberg, Berlin, New York, (2001). 

D. Braun and J. Martin, Nat. Commun. 2, 223 (2011). 

G. Lindblad, Commun. Math. Phys. 48, 119 (1976); V. 
Gorini, A. Kossakowski, E. C. G. Sudarshan, J. Math. 
Phys. 17, 821 (1976). 

O. Oreshkov and J. Calsamiglia, Phys. Rev. Lett. 105, 
050503 (2010). 

S. Nakajima, Prog. Theor. Phys. 20, 948 (1958); R. 
Zwanzig, J. Chem. Phys. 33, 1338 (1960). 

B. L. Hu, J.P. Paz, Y.H. Zhang, Phys. Rev. D 45 2843 

(1992). 

R. Karrlein, H. Grabert, Phys. Rev. E 55 153 (1997). 
F. Haake, R. Reibold, Phys. Rev. A 32 2462 (1985). 
J. H. An, and W. M. Zhang, Phys. Rev. A 76, 042127 

(2007) ; J. H. An, M. Feng, and W. M. Zhang, Quantum 
Inf. Comput. 9, 0317 (2009). 

C. -H. Chou, T. Yu, B. L. Hu, Phys. Rev. E 77 011112 

(2008) . 

J. P. Paz, A. J. Roncaglia, Phys. Rev. Lett. 100 220401 

(2008); Phys. Rev. A 79 032102 (2009). 

M. W. Y. Tu and W. M. Zhang, Phys. Rev. B 78, 235311 

(2008). 

M. W. Y. Tu, M. T. Lee, and W. M. Zhang, Quant. Info. 

Proc. 8, 631 (2009). 

J. S. Jin, M. T. W. Tu, W. M. Zhang, and Y. J. Yan, 
New J. Phys. 12, 083013 (2010). 

H. T. Tan and W. M. Zhang, Phys. Rev. A 83, 032106 
(2011). 

C. U Lei and W. M. Zhang, Ann. Phys. 327, 1408 (2012). 
J. Schwinger, J. Math. Phys. 2, 407 (1961); L. V. 
Keldysh, Sov. Phys. JETP, 20, 1018 (1965). 
H. J. Carmichael, An Open Systems Approach to 
Quantum Optics, Lecture Notes in Physics, Vol. ml8 
(Springer- Verlag, Berhn, 1993). 

H. P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, New York, 



11 



2007). 

[67] U. Weiss, Quantum Dissipativo Systems, 3rd ed. (World 
Scientific, Singapore, 2008). 

[68] L. P. KadanofF, G. Baym, Quantum Statistical Mechan- 
ics, (Benjamin, New York, 1962). 

[69] G. D. Mahan, many-Particle Physics, 3rd Ed. (Kluwer 
Academic/Plenum Publishers, New Yoek, 2000). 

[70] D. Chrusciriski and A. Kossakowski, Phys. Rev. Lett. 
104, 070406 (2010). 

[71] J. Piilo, K. Harkocn, S. Maniscalco, and K. -A. Suomi- 
nen, Phys. Rev. A 79, 062112 (2009). 

[72] H. N. Xiong, W. M. Zhang, X. G. Wang and M. H. Wu, 



Phys. Rev. A 82, 012105 (2010). 
[73] M. H. Wu, C. U. Lei, W. M Zhang, H. N. Xiong, Opt. 

Express. 18, 18407 (2010). 
[74] C. U Lei and W. M. Zhang, Phys. Rev. A 84, 052116 

(2011). 

[75] J. S. Jin, X. Zheng, and Y. J. Yan, J. Chem. Phys. 128, 
234703 (2008). 

[76] J. Maciejko, J. Wang, and H. Guo, Phys. Rev. B 74, 
085324 (2006). 

[77] T. Hayashi, T. Fujisawa, H. D. Cheong, Y. H. Jcong, and 
Y. Hirayama, Phys. Rev. Lett. 91, 226804 (2003). 



